Spin-lattice interaction in the quasi one-dimensional helimagnet LiCu2 02 
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The field dependence of the electron spin resonance in a helimagnet, LiCu202, was investigated 
for the first time. In the paramagnetic state a broad resonance line was observed corresponding to a 
g-factor of 2.3. In the critical regime around the paramagnetic to helimagnetic phase transition the 
resonance broadens and shifts to higher frequencies. A narrow signal is recovered at low temperature, 
corresponding to a spin gap of 1.5 meV in zero field. A comprehensive model of the magnons is 
presented, using exchange parameters from neutron scattering (T. Masuda et al. Phys. Rev. B 72 
014405), and the spin anisotropy determined here. The role of the quantum fluctuations is discussed. 

PACS numbers: 75.30.Ds, 75.25.+Z, 75.50.-y 
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Ground state solutions of the antiferromagnetic Hei- 
senberg model range from the Neel state to valence bond 
solids and resonating valence bonds, depending on the 
interactions, the dimensionality and connectivity of the 
system and other factors* 1 ^ When the ground state has 
long range Neel order, the excitation spectrum is gap- 
less at q = 0. This follows from the rotational invari- 
ance of the Hamiltonian: the creation of long wavelength 
excitations, where the relative angle of the neighboring 
spins changes only by an infinitesimal amount, costs very 
little energy. Gapless excitations, on the other hand, 
cause strong quantum fluctuations, suppressing or de- 
stroying the Neel stated The spin-orbit coupling leads to 
new terms in the Hamiltonian, including the single ion 
anisotropy, the exchange anisotropy, the Dzyaloshinskii- 
Moriya coupling and others. The coupling between the 
direction of the spins and the lattice removes the rota- 
tional invariance in many real materials, creates a gap in 
the spin wave spectrum, and contributes to the stability 
of the quasi-classical ground state. 

In this work we investigated the magnetic excitations 
in a quasi one-dimensional helimagnet, LiCu2C>2, in fields 
of T-14 T and at temperatures 2.5 K-60 K. The main 
result is that in the ground state two of the three magnon 
branches are gapped with A = 11.8 cm -1 =1.5 meV. Our 
calculations reveal that this gap is due to a spin-lattice 
coupling of D — 0.079 meV. We show that the proper 
treatment of the helical order of the spins explains the 
order-of-magnitude difference between the coupling pa- 
rameter and the gap, enabling this weak coupling to re- 
duce the quantum fluctuations significantly. 

LiCu2 02 crystallizes in an orthorhombic structure 
with the space group Pnma and lattice constants a = 
5.730 A, b = 2.86 A, c = 12.417 A 4 . The Cu 2+ ions 
carry magnetic moments with spin 1/2, and they form 
quasi one-dimensional chains with a zig-zag "ladders" 
along the b direction (Fig. 1). The material has two 
characteristics relevant to frustrated quantum magnets: 
the exchange interaction goes beyond the first neighbors, 
and the structure has a triangular motif susceptible to 




FIG. 1: (Color online.) Crystal structure of LiCu202, based 
on Ref. [j. The magnetic Cu 2+ ions are dark blue, the non- 
magnetic copper is green, the lithium is light blue and the 
oxygen is red. The dark blue bonds emphasize the triangular 
spin ladder. 



frustration. 

Quantum effects and fluctuations are expected to be 
important since the spins are 1/2 and the system is quasi 
one-dimensional. Indeed, in an early work Zvyagin et al. 
interpreted the broadening and disappearance of the ESR 
signal in terms of a "dimer liquid" state below the 24 K 
phase transition^. However, instead of a quantum liquid 
state, NMR, /iSK and neutron scattering results revealed 
a quasi-classical helical spin order along the chains, with 
a wavevector Q incommensurate to the lattic o 6 i 8 i 9 . All 
spins were found to be parallel to the a — b plane. In the 
q = ±Q region the "acoustic" magnon branches were 
observed, but the energy resolution was not sufficiently 
high to exclude or measure a possible gapAS. ESR has 
the resolution to address this question, and the applica- 
tion of the external static magnetic field can be used to 
explore the spin dynamics further. Note that in ESR the 
magnetic excitation is at q = 0, but in the presence of a 
magnetic superlattice it couples to the modes at q = ±Q 
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FIG. 2: Temperature dependence of the ESR signal at 12T 
magnetic field. The measured transmissions, relative to the 
transmission in zero field, 25K, are shown. Colors indicate the 
temperature. In the lower panel the same data is assembled 
into a two-dimensional intensity map. 



as well. 

The measurements were performed on an oriented sin- 
gle crystal sample of dimensions 3 mmx6 mm x 0.5 mm. 
According to the Laue diffractogram, the c axis is along 
the shortest dimension of the sample; the a and b direc- 
tions are parallel to the edges of the slab. There is a 
twinning in the a — b plane, as typical of most LiCu2 02 
samples. Spin resonance was detected in the transmission 
of the far infrared light, measured at Stony Brook Uni- 
versity's high magnetic field/infrared facility at the U12 
IR beamline of the National Synchrotron Light Source. 
The light propagated parallel to the static magnetic field, 
and passed through the sample along c direction. The po- 
larization of the incident light was controlled and set to 
several directions within the a — b plane. 

Figure 2 shows the temperature dependence of the spin 
resonance at 12 T field. The raw transmission curves 
were normalized to the transmission of the sample in T 
at 25 K, when no spin resonant absorption is expected. 
The oscillations seen in the baseline of the frequency de- 
pendence are residuals of the interference fringes seen in 
the raw spectra. These fringes are common and well un- 
derstood for samples of plane-parallel geometry. 

At high temperatures a broad resonance at the "free 
spin" frequency of hw — g/isH with a g factor of 



g = 2.3 is observed. The ESR line starts to broaden 
below 30 K; all of this is in agreement with earlier ESR 
investigationsiiii*i£. In the critical regime around 22 K- 
24 K the ESR line is about 3 T-4 T broad (full width at 
half maximum), and it shifts to higher frequency. Below 
15 K a narrower signal is recovered at a markedly higher 
frequency. This signal, characteristic to the development 
of a spin gap in the ordered phase, has not been seen in 
earlier spin resonance works£*ii*i£, due to limitations in 
accessible fields and frequencies in those studies. Mea- 
surements at several other fields (8 T,10 T and 14 T), 
yielded results qualitatively similar to the behavior pre- 
sented in Fig. 2, and will be published elsewhere. Choi 
et al. suggested a transition from helimagnetic to Neel 
orderi£ at T=9 K, but there is no indication of a change 
in the ESR signal in that range. 

Zvyagin and co-workers measured ESR in the param- 
agnetic state at 277 GHz (9.2 cm _1 )£. They observed 
the broadening and disappearance of the ESR signal, in- 
terpreting it in terms of the formation of a spin singlet 
liquid state above the ordering temperature. In fact, as 
our data shows, the integrated intensity of the signal re- 
mains approximately constant as we cross the phase tran- 
sition temperature. Vorotinov et al. observed an "an- 
tiferromagnetic resonance" (AFMR) signal at ^30 GHz 
(1 cm -1 )iAii2 in zero field, with an unexpected and unex- 
plained magnetic field dependence and anisotropy. There 
is no information about the intensity of this signal rel- 
ative to the ESR in the paramagnetic state. However, 
since nearly all of the spectral weight of the paramag- 
netic resonance remains in the high field/frequency range 
studied here, the AFMR signal at 1cm ~ 1 must be rela- 
tively weak, and may be related to low energy processes 
that are secondary to the physics of this material. 

Figure 3 shows the field dependence of the spin reso- 
nance at low temperatures. The three data sets in the up- 
per panel correspond to the different polarization states 
of the incident light. Each set of curves were obtained by 
dividing the measured spectrum with a reference spec- 
trum recorded at a temperature above the phase tran- 
sition in zero external field. The resonance lineshapes, 
frequencies and the intensities are similar for all three 
polarization states. The most striking feature is the dra- 
matic increase of the resonant absorption at higher mag- 
netic fields. 

To extract the field dependence of the resonant fre- 
quency we took the average of the three sets of data, and 
assembled an intensity map (Fig. 4). In terms of sim- 
ple, two-sublattice antiferromagnets both the intensity 
and the field dependence are surprising. The intensity of 
the antiferromagnetic resonance should have a very weak 
field dependence^ only. More importantly, the typical 
field dependence of the resonance frequency in antiferro- 
magnets (u> oc y/ ' H 2 + Hq) can not be fitted to the data 
without residual systematic deviations. 

The model used for the magnon spectrum^ was 
adopted to evaluate the results presented in Fig. 3. We 
extended it to finite magnetic fields perpendicular to the 
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plane of the spins, and added the terms responsible for 
the spin gap. The Hamiltonian is: 

TL = JiSijSi+i.j + J2SijSi+2,j + JiS^jSi+i.j + 

i,3 

+J±SijSi, j+1 - g^BHSfj + TL' , (1) 

where 

H' = D s {Sl 3 f or W ILS'!^!..., (2) 

The indices i and j run along and perpendicular to the 
double chains, respectively. J\ is the coupling between 
along the diagonal "rungs" of the ladder, J 2 and J 4 
are the nearest neighbor and the second neighbor cou- 
plings along the chain, and J± is the interchain coupling. 
The "easy plane" is represented either by the single ion 
anisotropy of D s > 0, or by the exchange anisotropy of 
D ex < (the negative D ex is required because the ori- 
entation of the nearest neighbor spins has an antiferro- 
magnetic character). Following conventions established 
in the literatur e) 16 ! 17 , the reference frame for the spin 
components is selected so that the magnetic field H is 
applied along the y directions and in zero field the spins 
are in the x — z plane (this makes the crystallographic c 
direction parallel to the y axis). In finite magnetic field 
the spin directions are on a cone whose axis is parallel to 
V- 

The ground state and the spin wave excitations of 
this m odel were found with the methods described in 
Refs. Ildl6ll7l Following Masuda et alH we used 
Ji = 6.4 meV, J 2 = -11.9 meV, J 4 = 7.6 meV and 
J_L = 1.8 meV, and in zero field and with no anisotropy 
we reproduced the fits to the published spin wave spec- 
trum. After carefully tracing the factor 2's, we believe 
that the coupling constants published in the preprint ver- 
sion of Ref. [l0| were correct. 

In general, a helical spin arrangement possesses 3 Gold- 
stone modes, one at w(0) corresponding to the free ro- 
tation of the ordered moments within the helical plane, 
and two degenerate ones at w(±Q) associated with the 
tilting of the plane of the helix 17,18 . The "easy plane" 
anisotropy generates a finite spin gap at ±Q, where Q = 
(7r/a, </>/&) is the ordering wave vector with <j> = 5.2rad 
(the angle between subsequent spins on the same leg of 
the ladder). In terms of the Fourier transform of the ex- 
change interaction, J(q) = 2( J\ cos(6gb/2) + J 2 cos(fr<2&) + 
J 4 cos(26qh) -I- J± cos(aq a )), the gap is 

A = 2Sy/j'D s or A = 2S\/ J' cos(0/2)L> ex (3) 

for single ion or exchange anisotropy, respectively. The 
effective exchange coupling J' is defined as 

= J(2Q) + J(0) _ J(Q) 
'4 2 [ 1 

The application of external magnetic field lifts the de- 
generacy of the ±Q points, resulting in distinct spin gaps 
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FIG. 3: Field dependence of the transmission at 2.5K. The 
three sets of curves represent the transmission relative to the 
transmission at high temperature and zero field. Top: par- 
tially polarized light, relative to 25K. Middle: light polarized 
parallel to the long axis of the sample, relative to 20K. Bot- 
tom: polarization along the short axis, relative to 25K. 
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FIG. 4: Upper panel: Intensity map based on the average of 
the three datasets shown in Fig. 3. Lower panel: Intensity 
map generated from the calculated frequencies and suscepti- 
bilities, based on the model described in the text. The lines 
represent the resonance frequency corresponding to the two 
branches, also plotted in the upper panel. 



and different spin susceptibilities (Fig. 4). The contin- 
uous lines in the Figure represent the center-frequencies 
of the two branches generated by the application of the 
external field. The strength of the resonance line was 
set according to the calculated spin susceptibility, and 
lifetime effects were simulated by Lorentzian lineshapes 
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of O.lmeV relaxation rate, independent of field. On the 
lower branch, the strong absorption at high fields grad- 
ually decreases as the field is lowered, due to the field 
dependence of the susceptibility. The susceptibility be- 
longing to the upper branch is so small that the the cor- 
responding resonance signal is below the noise level. 

In order to obtain the near-perfect agreement with the 
experiment we adjusted two parameters: the magnitude 
of the anisotropy, D s = D ex cos(0/2) = 0.079, was set 
to reproduce the zero-field gap, and the <?-factor was 
g = 1.85. Instead of using a single ion anisotropy in Eq. 
l[TJl , one may use exchange anisotropy or Dzyaloshinskii- 
Moriya interaction, with very similar results^. The mi- 
croscopic origin of all of these interactions is the spin 
orbit coupling, but ab-initio calculations are difficult. 
Moriya 19 estimated the anisotropic component of the ex- 
change interaction in the order of D ~ (Ag/g) 2 J, where 
Ag is the g-factor shift in the paramagnetic state and 
Ag/g is a measure of the spin orbit coupling. A spin- 
lattice coupling in the range of D ~ 100 is there- 
fore consistent with the g-factor ranging between 1.95 
and 2.3 in the paramagnetic stateiii*i£ The energy gap 
scales linearly with the spin orbit and exchange cou- 
plings, (A ~ Ag/g) J. 

Quantum fluctuations suppress the sublattice mag- 
netic order by the amount of AS* = J ^pjd D q (here D is 
the dimensionality of the lattice, w(q) is the magnon ex- 
citation energy at wavevector q, and g(q) is expressed in 
terms of J(q)) Ail. With no spin-lattice coupling we ob- 
tained AS* = 0.150, a 30% drop relative to spin S = 1/2, 
which is to be contrasted with A5* = 0.197 for a two di- 
mensional antiferromagnet. In the presence of the easy 
plane term the corresponding value is AS* = 0.125. The 
coupling to the lattice makes the quasi-classical state 
more stable. Further reduction of AS is seen in the 
presence of static magnetic field. Further evidence for 



reduced sublattice magnetization follows from g-factor 
determined in the fit. The macroscopic energy scale of 
the interaction with the magnetic field is determined by 
the product gfj,BHS e ff. In our treatment the quantum 
effects were neglected, and instead of a reduced S the 
quantum correction appears as a reduction of g. 

In conclusion, we established that the spin wave spec- 
trum of LiCu202 has a gap of 1.5meV, and interpreted 
the experiments in terms of a model including an easy 
plane anisotropy. This term in the Hamiltonian accounts 
for the experimentally observed spin direction in the crys- 
tallographic a — b plane, and it also contributes to the 
stability of the quasi-classical ground state. There are 
intriguing possibilities for similar measurements on heli- 
magnets with external field applied in the plane of the 
spin rotation. The theoretical methods applied in this 
work can not be readily adapted to this new configura- 
tion. It is known that the static spin order turns into a 
a soliton latticej2L2£ but calculations for the ESR exci- 
tations have not been performed yet. 
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